home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DGECO.z / DGECO
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDGGGGEEEECCCCOOOO((((3333FFFF))))                                                            DDDDGGGGEEEECCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DGECO   - DGECO factors a double precision matrix by Gaussian elimination
  10.      and estimates the condition of the matrix.
  11.  
  12.      If  RCOND  is not needed, DGEFA is slightly faster.  To solve  A*X = B ,
  13.      follow DGECO by DGESL.  To compute  INVERSE(A)*C , follow DGECO by DGESL.
  14.      To compute  DETERMINANT(A) , follow DGECO by DGEDI.  To compute
  15.      INVERSE(A) , follow DGECO by DGEDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE DGECO(A,LDA,N,IPVT,RCOND,Z)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAA DOUBLE PRECISION(LDA, N)
  25.         the matrix to be factored.
  26.  
  27.      LLLLDDDDAAAA INTEGER
  28.         the leading dimension of the array  A .
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  On Return
  32.  
  33.      AAAA an upper triangular matrix and the multipliers
  34.         which were used to obtain it.
  35.         The factorization can be written  A = L*U  where
  36.         L  is a product of permutation and unit lower
  37.         triangular matrices and  U  is upper triangular.
  38.  
  39.      IIIIPPPPVVVVTTTT INTEGER(N)
  40.         an INTEGER vector of pivot indices.
  41.  
  42.      RRRRCCCCOOOONNNNDDDD DOUBLE PRECISION
  43.         an estimate of the reciprocal condition of  A .
  44.         For the system  A*X = B , relative perturbations
  45.         in  A  and  B  of size  EPSILON  may cause
  46.         relative perturbations in  X  of size  EPSILON/RCOND .
  47.         If  RCOND  is so small that the logical expression
  48.         1.0 + RCOND .EQ. 1.0
  49.         is true, then  A  may be singular to working
  50.         precision.  In particular,  RCOND  is zero  if
  51.         exact singularity is detected or the estimate
  52.         underflows.
  53.  
  54.      ZZZZ DOUBLE PRECISION(N)
  55.         a work vector whose contents are usually unimportant.
  56.         If  A  is close to a singular matrix, then  Z  is
  57.         an approximate null vector in the sense that
  58.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  59.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDGGGGEEEECCCCOOOO((((3333FFFF))))                                                            DDDDGGGGEEEECCCCOOOO((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      Subroutines and Functions LINPACK DGEFA BLAS DAXPY,DDOT,DSCAL,DASUM
  75.      Fortran DABS,DMAX1,DSIGN
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.